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ABSTRACT 


The  experimental  results  for  convective  heat  transfer  and 
flow  friction  characteristics  of  some  screen,  plate-fin  and 
skewed  fin  matrices  of  various  geometry  are  presented.  The  heat 
transfer  data  were  obtained  using  the  transient  technique. 

The  results,  which  are  presented  in  graphic  and  tabular 
form,  should  be  of  practical  use  for  compact  heat  exchanger 
design. 

The  results  of  a  digital  computer  solution  giving  the 
maximum  slope  of  the  generalized  heating  curve  as  a  unique 
function  of  NTU  are  presented  in  tabular  form  in  Appendix  II, 
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INTRODUCTION 


There  are  two  courses  open  to  the  designer  of  a  compact 
heat  exchanger  for  the  transfer  of  a  fixed  heat  load  consistent 
with  an  allowable  pressure  drop  utilizing  a  new  type  of  matrix 
arrangement  that  is  not  covered  in  present  heat  transfer  litera¬ 
ture.  A  matrix  will  be  defined  as  a  porous  solid  possessing  a 
high  ratio  of  void  surface  area  to  bulk  volume,  and  a  void  geo¬ 
metry  which  permits  the  flow  of  gas  or  a  liquid  through  it. 

The  first  course,  a  trial  and  erx'or  technique,  consists  of 
analytically  predicting  the  heat  transfer  and  friction  character¬ 
istics  of  the  matrix.,  designing  and  testing  a  scale  version  of  the 
heat  exchanger,  and  modifying  the  design  until  the  required  perfor¬ 
mance  is  achieved.  The  time  and  cost  involved  in  utilizing  this 
method  is  generally  unacceptable  to  the  designer. 

The  second  method  consists  of  generating  the  required  heat 
transfer  and  friction  data  on  a  small  module,  and  utilizing  the 
results  for  the  full  scale  design.  This  method  is  quite  flexible, 
since  many  different  matrix  arrangements  can  be  tested  at  mini¬ 
mum  cost,  and  an  optimum  surface  chosen. 

One  of  the  satisfactory  methods  of  measuring  the  convective 
heat  transfer  characteristics  of  porous  solids  utilizes  the  so- 
called  transient  test  technique.  This  consists  of  introducing  a 
step  function  in  the  temperature  of  the  fluid  passing  through  the 
matrix,  which  is  initially  at  uniform  temperature,  and  recording 
the  matrix  temperature-time  history.  Using  the  results  of  the 


analysis  developed  by  d.  L»  Locke  (l)’*’  for  the  determination 
of  NTU  from  the  generalized  heating  curve  of  maximan  tempera¬ 
ture  gradient  versus  time,  it  is  then  possible  to  compute  the 
Stanton  number.  The  'Fanning'  friction  factor  is  obtained  from 
the  pressure  drop  across  the  matrix. 

The  objectives  of  this  thesis  were: 

a.  to  design,  fabricate  and  evaluate  the  performance 
of  a  test  section  to  utilize  the  transient  method  fpr  deter¬ 
mining  heat  transfer  coefficients  of  small  modules,  and 

b,  to  obtain  the  heat  transfer  coefficients  and  fluid 
frictional  characteristics  of  screen,  triangular  fin  and  plate- 
fin  matrices  of  various  geometry. 


^Numbers  in  parentheses  refer  to  Bibliography  on  page  21. 
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EXPEK  MENTAL  APPARATUS 


A  schematic  diagram  of  the  apparatus  is  shown  in  Fig,  8. 
Basically,  the  system  consists  of  a  plastic  section  with  a 
sliding  drawer  for  holding  the  test  matrix,  a  flow  measuring 
system,  a  pressure  measuring  system,  a  compressor  to  provide  flow 
through  the  test  section,  an  air  heater,  a  temperature  measuring 
system  and  piping. 

Flow  Rate  Control 

The  compressor  is  run  at  a  constant  rate  of  approxi¬ 
mately  300  c.f.m.  A  sliding-gate  valve  on  a  T-section  of  the 
inlet  piping  of  the  compressor  is  used  for  coarse  control  of 
air  flow  through  the  test  section.  There  is  also  a  fine  control 
valve  for  making  minor  adjustments  to  the  air  flow  rate.  A  thin 
plate  orifice  section  for  flow  measurement,  designed  according  to 
A.S.M.E.  specifications  (3),  is  included  in  the  inlet  piping. 
Hatrix  Heating  System 

A  small  blower,  with  a  heating  element  controlled  by  a 
115 -volt,  60-cycle  AC  'Variac',  is  positioned  to  heat  the  matrix 
when  the  drawer  is  removed  from  the  plastic  test  section.  A 
transitional  section,  with  flow  stabilizing  screens,  is  attached 
to  the  blower  to  accomplish  uniform  heating  of  the  matrix. 
Temperature  Measuring  System 

A  five-thei'mocouple  upstream-temperature  probe  is  located 
in  the  center  of  the  pipe  approximately  midway  between  the  ori¬ 
fice  and  test  section.  The  upstream-downstream  air  temperature 
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difference  is  measured  by  these  thermocouples,  connected  in 
series  so  as  to  oppose  the  3  center  thermocouples  of  a  13= 
thermocouple  tray  located  immediately  after  the  matrix. 

The  thirteen-thermocouple  tray  may  be  used  indepen¬ 
dently  to  find  the  temperature  distribution  of  the  downstream 
air.  all  thermocouples  used  are  30-gauge  irozL-constantan, 

The  temperature  difference  and  distribution  is  recorded  on  a 
Leeds  &  Northrup  'Speedomax',  Model  G,  emf  recorder  with  a 
variable  range  of  from  1-20  inv. 

Ambient  air  temperature  is  read  on  a  mercury  thermo¬ 
meter  located  at  the  inlet  of  the  entrance  piping. 

Sliding  Drawer  Matrix  Holder 

The  sliding  drawer  matrix  holder  is  constructed  of 
half-inch  thick  plastic.  It  is  lined  with  balsa  wood,  al¬ 
lowing  for  a  3o08-inch  diameter  passage  for  holding  the  screen 
matrix,  when  testing  the  triangular  fin  matrices,  a  3«I23~ 
inch  square  passage  is  provided  for  holding  the  matrix. 
pressure  Measurement  System 

The  pressure  drop  across  the  orifice  is  measured  on  a 
lJj-inch  U-tube  manometer  filled  with  water.  The  pressure  drop 
across  the  matrix  is  measured  on  either  a  30-inch  differential 
manometer  or  an  inclined  160  mm.  manometer  filled  with  water. 

Ambient  air  is  drawn  through  a  3.08-inch  diameter  pipe 
to  the  test  section.  When  testing  the  triangular  fin  matrices. 


h 


a  transition  piece  to  a  3 •125-inch  square  was  inserted  before  and 
after  the  test  section,,  Galvanized  ducting  is  used  between  the 
test  section  outlet  and  compressor  inlet. 
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METHOD 


A.  Summary  of  Theory 

The  theoretical  analysis  of  the  transient  heat  transfer 
behavior  of  a  porous  medium  was  first  presented  by  3 human n  in 
1929.  Shumann  considered  a  homogeneous  porous  medium  at  a  uni- 
foxin  temperature  through  which  a  fluid  of  the  sane  temperature 
is  flowing.  At  a  certain  instant,  the  temperature  of  the  enter- 
ing  fluid  is  assumed  to  change  to  a  higher  or  lower  value.  The 
problem  is  then  to  find  the  temperature  of  the  fluid  and  solid 
as  functions  of  time  and  position  in  the  matrix. 

The  following  idealizations  and  initial  and  boundary 
conditions  have  to  be  satisfied  in  order  to  simplify  the 
problem  2 
Idealizations 

1.  The  fluid  specific  heat  and  viscosity  are  constant. 

2.  The  thermal  conductivity  of  both  fluid  and  solid  is 
zero  in  the  direction  of  flow  and  is  infinite  within 
the  fluid  and  solid  in  the  direction  nonrial  to  the  flow. 

3.  The  flow  is  steady. 

li.  The  porous  solid  is  homogenous. 

Initial  and  Boundary  Conditions 

1.  At  time  9=0,  the  temperature  of  the  entering  fluid 
changes  instantaneously. 

2.  Initially,  the  core  is  at  a  uniform  temperature. 

3.  Ho  heat  passes  the  core  boundaries. 
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Shurnann  found  the  solutions  for  the  fluid  and  metal 
temperatures  to  be  2 


and 


Shurnann' s  solution  requires  that  a  theoretical  curve  be 
generated  which  will  match  the  experimental  curve  to  deter¬ 
mine  the  heat  transfer  coefficient. 

In  1950,  Locke  (1)  developed  a  method  in  which  only 
the  maximum  slope  of  the  theoretical  and  experimental  curves 
need  comparison.  This  method  eliminates  the  requirement  of 
a  large  number  of  theoretical  carves  and  the  possibility  of 
error  due  to  the  displacement  of  the  experimental  heating 
curve  either  vertically  or  horizontally.  Locke  derived  an 
expression  for  the  slope  of  the  generalized  heating  curve 
for  the  fluid,  (equation  1  above),  finding  that  at  X'-L  for 
which  2  becomes  equal  to  the  MTU: 


From  this,  it  was  found  that  the  maximum  slope  of  the  generalized 
heating  curve  is  a  unique  function  of  MTU.  Thus,  only  the  maximum 
slope  of  the  experimental  heating  curve  need  be  found  in  order  to 
obtain  the  NTU  and,  thereby,  the  Stanton  number.  This  method  has 
been  used  by  several  investigators,  and  their  results  have  been 
compiled  by  Kays  &  London  (2). 

■a- 

Locke,  in  his  solution  of  Equation  3,  resorted  to  an 
approximation  of  the  modified  Bessel  function  which  in  the  low 
range  of  the  argument  was  subject  to  large  error.  A  solution 
to  Equation  3  has  been  carried  out  by  Lt.  J.  A.  Inglis,  RON, 
using  the  exact  infinite  series  of  the  modified  Bessel  function, 
on  a  CDC  160U  digital  computer.  The  results  of  this  solution 
are  given  in  Appendix  II  and  were  the  ones  used  by  the  author 
in  reducing  his  experimental  data. 

B.  Experimental  Techhique 

Ambient  air  is  drawn  through  the  inlet  pipe  and  allowed 
to  flow  through  the  test  matrix,  while  the  flow  rate  is  ad» 
justed,  after  constant  flow  is  established,  the  sliding  drawer 
containing  the  matrix  is  pulled  out  of  the  test  section  and  into 
the  heated  air  stream. 

When  the  matrix  temperature  is  steady,  as  indicated  by 
monitoring  thv.  thermocouples,  the  drawer  is  rapidly  inserted 
into  the  ambient  air  stream  and  the  expei  iment  is  started. 

The  ambient  air  temperature,  orifice  pressure  drop  and 
matrix  pressure  drop  are  recorded,  and  the  transient  upstream 
to  downstream  temperature  difference  is  recorded  on  the  Leeds 
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&  Northrup  'Speedomax1  Recorder,  After  the  temperature  differ- 
ence  has  returned  to  zero,  the  run  is  terminated „ 

The  following  data  are  recorded: 

a,  atmospheric  pressure 

b,  atmospheric  temperature 

c.  upstream  to  downstream  temperature  difference 

d.  pressure  drop  across  the  orifice 

e.  pressure  drop  across  the  matrix 

f,  orifice  diameter 

In  this  investigation,  the  temperature  difference  between 
the  inlet  fluid  and  the  downstream  fluid  temperature  was  recorded 
directly  so  the  maximum  slope  could  be  determined  from  the  ex¬ 
perimental  trace  with  no  replotting.  For  a  complete  set  of 
sample  calculations,  see  Appendix  I. 

C.  Presentation  of  Results 

To  avoid  the  confusion  often  encountered  when  a  large 
number  of  arbitrarily  defined  parameters  is  used,  the  same 
definitions  will  be  applied  throughout  to  the  various  matrix 
geometries. 

The  basic  data  and  results  for  each  surface  are  pre¬ 
sented  in  both  tabular  (Tables  I  -  VI)  and  graphic (Figs.  1-5) 
forms,  as: 

KStKPr^3  =  *2  (“Re) 

f  =  62(NRe) 
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The  Reynolds  number  is  based  on  a  hydraulic  diameter, 


defined  as  follows: 

Ur  G 
NKe=  4  h 


with 

^rh  ■  Dh  -  UAcn/A,  ft., 

G  =  Wf/Ac,  IV m  hr-ft  2 

The  Stanton  number  is  defined  as: 


NSt  = 


_  h 


GCf 


-  NTII  A_/A 


where 


NTU  =  ^/\i£C£  -  number  of  heat  transfer  units  of  the 
matrix. 

The  friction  factor  is  defined  on  tine  basis  of  an  equivalent 
shear  force  in  the  flow  direction  per  unit  of  heat  transfer  area. 
Entrance  and  exit  losses  are  included  in  the  friction  factor, 
f  =  ^  matrix  2^c  Ac/C:  v.fA 

Actually,  Prandtl  number  was  not  a  test  variable,  as  the 
working  fluid  was  in  all  cases  air  at  moderate  temperatures. 

Ove'  a  moderate  range  of  Prandtl  numbers,  the  effect  of  Prandtl 
number  to  the  2/3  power  is  used  as  an  approximation. 


D.  Limitations  and  Uncertainties 

Gare  must  be  used  to  approach  the  idealizations  and  initial 
and  boundary  conditions  of  the  theoretical  analysis,  which  are 
restated  belovj: 
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Idealizations 


1„  The  fluid  specific  heat  and  viscosity  are  constant. 

2.  The  thermal  conductivity  of  both  fluid  and  solid  is 
zero  in  the  direction  of  flow  and  is  infinite  within 
the  fluid  and  solid  in  the  direction  normal  to  the  flow. 

3 »  The  flow  is  steady. 

h.  The  porous  solid  is  homogenous . 

Since  the  fluid  used  is  air,  the  first  idealization  can  be 
approached  if  temperature  differences  are  small.  As  the  maximum 
temperature  difference  is  approximately  20°  F. ,  the  arithmetic 
average  will  be  plus  or  minus  1,0$  in  viscosity  and  negligible  in 
specific  heat. 

As  for  the  second  idealization,  it  is  recognized  that  ther¬ 
mal  conduction  along  the  flow  length  occurs  and  increases  as  the 
flow  rate  is  lowered.  Its  effect  can  be  lowered  by  shortening, 
or  dividing  the  length  of  matrix  in  the  flow  direction.  Mondt  (?) 
indicates  that  as  the  Reynolds  number  is  decreased  the  conduction 
effect  will  be  small  below  an  NT'J  of  11.79,  have  no  effect  at 
NTU  =  11.79,  and  increase  appreciably  as  the  Reynolds  number  is 
further  reduced. 

In  the  direction  normal  to  the  flow,  the  conduction  is 
negligible  since  the  temperature  gradient  at  any  section  is 
essentially  zero.  This  gives  a  condition  which  is  equivalent  to 
infinite  conductivity. 

The  third  idealization  is  approached  by  running  the  com- 
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pressor  at  rated  capacity  for  all  runs,  so  that  pulsing  does  not 
occur. 

The  fourth  idealization  is  approached  by  careful  selection 
of  the  porous  solid  to  be  tested. 
initial  and  Boundary  Conditions 

1.  at  time  9=0,  the  temperature  of  the  fluid  flowing 
through  the  matrix  changes  instantaneously. 

2.  Initially,  the  core  is  at  a  uniform  temperature. 

3 o  No  heat  passes  the  core  boundaries. 

The  first  condition  can  be  approached  by  rapidly  intro¬ 
ducing  the  heated  matrix' into  the  flow  stream  of  constant- 
temperature  ambient  air. 

The  second  condition  is  satisfied  by  placing  screens  in 
the  heater  duct  to  produce  a  uniform  velocity  profile  of  the 
heated  air.  The  downstream  side  of  the  matrix  Ls  monitored 
by  a  thirteen- thermocouple  tray  to  check  for  a  uniform  tem¬ 
perature  distribution  over  the  section. 

The  third  condition  can  only  be  minimized  by  careful 
design  of  the  matrix  holder.  To  lower  the  transfer  of  heat 
from  the  matrix  holder  to  ti  e  matrix  on  cooling,  the  matrix 
holder  is  lined  with  balsa  wood . 

The  uncertainties  of  physical  constants  in  the  matrix 
and  duct  measurements  are  listed  below: 
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Quantity 

cS 

NPr 

cf 

W„ 


Uncertainty  ai  20:1  odds 

“U  ,o% 

±1,2% 

-0.h% 

-1.0/1 

negligible 


toA% 


A 


-l.C# 

-l.O^I 


Instrumentation  uncertainties  follow: 


Quantity 

ti-tfl 


Uncertainty  at  20:1  odds 

“0.5:1 


^orifice  -1.0$ 

^matrix  ^2.0^ 

Following  is  a  summary  of  the  uncertainty  analysis  of 
the  Reynolds  number,  friction  factor,  NTU  and  Stanton  number 
For  Repeatability' 


Quantity 

%e 

f 

NTU 


N 


St 


N3tNpr" 


/3 


Uncertainty  at  20:1  odds 

*2.1# 

-3.3  < 

i3o0-ior' 

is.o-io: 

-8.0-10,,. 
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The  NTU,  and  N3^Npr2/3  are  represented  with  a 
variable  range  of  uncertainties  because  NTU  does  not  vsxy 
linearly  with  the  maximum  slope  of  the  cpoling  curve. 


NTU 

Uncertainty  at 

6 

iicy 

10 

~  9.0% 

16 

8 . 0% 

20 

~  9.0% 

Reproducabi] ity 

Quantity 

Uncertainty  at 

Nfte 

-  3c 0f% 

f 

* 

NTU 

-13-22)6 

Not 

-18-22^ 

N3tNPr2/3 

-13-22$ 

NTU 

Uncertainty  at 

6 

-22$ 

10 

«20# 

16 

-18$ 

20 

-19% 

1U 


DJSSUSSMN  OF  RESULTS 


To  evaluate  the  performance  of  the  test  section, 

60  x  6C  x  0.0070  and  10  x  10  x  0.025  13/3  stainless  steel 
screen  matrices  which  .'/or :•  similar  to  those  used  by  other 
investigators  were  tested  (see  Table  711  and  Fig,  9).,  Fig.  1 
compares  the  results  of  the  present  investigation  with  the  be 
interpretation  of  the  data  obtained  by  Coppagc  and  Tong  (9) 
also  on  10  x  1C  x  0.025  screen  matrices.  Fig.  2  compares  the 
results  of  the  present  inv  ?stigation  of  60  x  60  x  C.0Q70 
screens  with  the  best  interpretation  of  the  data  obtained  by 
Cop  page  and  Tong  on  60  x  60  x  0,0075  screens.  Th'1  heat 
transfer  rnd  fluid  friction  characteristics  obtained  for  the 
screen  matrices  are  well  vithir.  the  range  of  uncertainty  of 
reproducability  of  Coppage,  Ton,  and  the  author. 

In  addition  to  the  screen  matrices,  which  were  tested 
primarily  to  evaluate  the  performance  of  the  test  section, 
three  different  geometry  matrices  were  investigated:  two  of 
the  plate-fin  type  (Fig.  10)  and  one  fin  type  (Fig.  11)  in 
which  the  elements  were  skewed  to  prevent  meshing  of  oho  fins 
Table  VIII  lists  the  details  of  these  matrices, 

Results  for  the  Type  U30  stainless  steel  plate-fin 
matrix  are  shown  in  Fig.  3  and  listed  in  Tables  III  and  IV. 
From  th*  figure  it  is  3een  that  in  the  higher  Reynolds  number 
range,  the  heat  transfer  behavior  follows  the  characteristic 
straight  line  with  a  clone  of  slightly  less  than  t.h*  minus  on 
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*  t,  a  Reynolds 


expected  tor  laminar  flow  in  smooth  channels, 
number  of  about  100  and  a  test  NTU  of  about  lj>,  the  data  begins 
to  depart  from  the  straight  l.ine  behavior,  As  the  Reynolds  num¬ 
ber  decreases,  this  departure  becomes  more  pronounced  and  finally 
reaches  a  point  where  it  a  i  pears  that  the  heat  transfer  coef¬ 
ficient  becomes  constant  or  independent  of  the  Reynolds  number. 

This  behavior  is  due  to  the  longitudinal  beat  conduction  within 
the  matrix,  and,  as  mentioned  earlier,  was  expected  to  start 
showing  up  when  the  MTU  theoretically  had  a  value  of  11.77(7). 

Since  the  conduction  parameter'”  for  a  given  geometry 
and  material  varies  inversely  as  the  Reynolds  number,  the 
smaller  the  Reynolds  number  the  greater  the  effect.  The  con¬ 
duction  parameter  also  varies  inversely  as  the  flow  length, so 
that  a  first  observation  -would  indicate  that  decreasing  th<° 
length  does  not  improve  the  situation;  but,  since  the  NTU  is 
directly  proportional  to  the  heat  transfer  area,  decreasing  the 
matrix  length  sufficiently  will  decrease  the  NTU  to  allow  test¬ 
ing  a  lower  Reynolds  number  range.  In  thin  investigation  the 
matrix  length  was  not  decreased  because  it  was  only  desired  to 
establish  the  straight  line  laminar  ■'low  behavior  and  sufficient 
data  was  obtained  to  do  this  without  altering  the  matrix  length. 
Raving  established  the  heat  transfer  behavior  in  the  higher  laminar 
Reynold.*?  number  range  thus  allows  extrapolation  of  the  data  into  the 
lower  Reynolds  number  range.  It  should  be  noted,  however,  that  when 

\  f 

,fnhfi  conduction  parameter  is  define* A  as  ~_s_ 


such  information  as  this  is  utilized  in  heat  exchanger  design, 
conduction  plays  an  important  part  in  the  heat  exchanger's 
performance,  and  there  has  been  no  attempt  in  this  investigation 
to  obtain  correlations  which  would  predict  the  effects  of  con- 
duction. 

The  friction  factor  curve  for  the  Type  U30  stainless 
steel  plate  fin  matrix  shows  essentially  a  straight  line 
behavior  with  a  slope  approximately  that  of  the  heat  transfer 
curve  as  would  be  expected  for  laminar  flow  conditions. 

Results  for  the  aluminum  plate-fin  matrix  are  shown 
in  Fig,  li  and  Table  V,  For  this  matrix  the  only  geometry  change 
was  that  of  the  hydraulic  diameter,  which  was  accomplished  by 
reducing  the  fin  height  and  pitch.  This  change,  along  with  the 
use  of  thinner  material,  provided  for  a  greatly  increased 
ratio  of  heat  transfer  area  to  unit  volume  and  of  free  flow  area 
to  frontal  area  compared  with  the  Type  ii30  stainless  steel  plate- 
fin  matrix.  The  experimental  results  show  the  same  type  behavior 
as  previously  discussed  with  the  conduction  effect  beginning  at 
a  Reynolds  number  of  about  200  where  the  test  NTU  was  about  13, 

In  Fig,  £  and  Table  VI  are  shown  the  results  for  the 
more  radical  of  the  geometry  changes.  For  this  matrix  the  plates 
were  removed  leaving  only  the  fins;  and  to  prevent  meshing  of 
the  fins,  the  elements  were  made  so  that  the  total  angle  between 
two  elements'  flow  channels  was  25°,  This  angle  was  chosen 
simply  because  it  was  an  easy  one  to  fabricate  and  sufficiently 
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removed  from  a  zero  angle  to  produce  a  large  effect.  It  was 
believed  that,  with  such  an  obviously  tortuous  flow  channel  the 
shape  of  the  result  curves  would  more  nearly  approach  those  of 
the  screen  matrices.  That  this  was  not  the  case  is  readily 
evident  by  comparing  Fig.  5  with  either  Fig.  1  or  2.  Instead, 
the  behavior  of  this  matrix  was  quite  similar  to  that  of  the 
plate-fin  type  and  will  be  discussed  in  more  detail  later.  The 
conduction  effect  is  seen  to  begin  at  a  Reynolds  number  of  about 
3C0  for  which  the  test  MTU  was  about  12. 

For  comparison  purposes  the  best  interpretation  of  the 
heat  transfer  characteristics  of  the  three  finned  matrices  have 
been  replotted  and  extended  in  Fig.  6,  and  the  fluid  friction 
characteristics  reproduced  in  Fig,  7. 

Comparing  the  stainless  steel  plate-fin  and  aluminum 
plate-fin  matrices  shows  the  effect  of  reducing  the  hydraulic 
diameter  by  reducing  size  and  pitch.  The  smaller  finned 
aluminum  matrix  has  approximately  the  same  fluid  friction 
characteristics  as  the  stainless  '••tee  l  matrix,  while  the  heat 
transfer  characteristics  are  raised  bl%>  This  would  indicate 
that  the  fins  in  plate-fin  heat  exchangers  could  be  reduced  in 
size  to  the  limit  allowed  by  considerations  other  than  heat 
transfer. 

The  comparison  of  the  aluminum  plate-fin  and  skewed  fin 
matrices  of  the  same  fin  size  shows  the  effect  of  removing  the 
plate  and  skewing  by  25°,  These  changes  increase  the  fluid 
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friction  characteristics  by  \00%  and  the  heat  transfer 
characteristics  by  73o^/»  Although  thn  heat  transfer  rate  is 
increased,  this  nay  be  an  expensive  means,  especially  in  gas 
turbini  regenerators  where  a  1,*  increase  in  pressure  drop  due 
to  flow  friction  costs  more  than  a  1 %  increase  in  specific  fuel 
consumption.  Possibly  there  is  a  skew  angle  at  which  the  in- 
crease  in  heat  transfer  characteristics  is  sufficiently  greater 
than  the  increase  in  friction  characteristics  to  warrant  its 
use. 

Comparing  the  stainless  steel  plate-fin  and  skewed 
aluminum  fin  matrices  shows  an  increase  of  \$9%  in  the  heat 
transfer  characteristics  and  an  increase  of  as  much  as  120^  in 
the  fluid  friction  characteristics.  In  this  case,  the  heat 
transfer  area  per  unit  volume  is  nearly  the  same  for  these 
surfaces,  but  the  hydraulic  diameter  for  the  skewed  fins  is 
practically  $0jo  greater,  indicating  that  there  would  be  less 
chance  of  fouling  with  the  skewed  arrangement. 


ic; 


JONJIUSajNS 


1.  A  successful  lest  section  has  been  designed  and  fabricated 
to  obtain  convective  heat  transfer  and  friction  coefficients 
for  small  modules. 

2.  Convective  heat  transfer  and  friction  coefficients  have  been 
obtained  for  some  screen,  triangular  fin  and  plate-fin 
matrices  of  various  geometries. 

3«  The  experimental  results  of  heat  transfer  behavior  and 
friction  factors  for  the  plate-fin  and  fin  matrices, 
summarized  in  Figs.  '6  and  7,  indicate  that  there  may  be 
better  geometries  for  compact  heat  exchanger  applications 
than  plate-fin  types  for  comparable  ratios  of  heat  transfer 
area  to  unit  volume. 
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.2/3  and  f^Versus  Npe  for  10x10  x  .025  18/8  stainless  eteel 
Figure  1 
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JStNPr2^3  and  f^  versus  Npd  for  60x60  X  .007  10/6  stainless  steel 
Figure  2 
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IJstNpr2/5  and  f  versus  NBe  for  aluminum  plata-fin 
r  Figure  4 


H  2/3  and  f  versus  SRo  for  25°  skewed  aluminum  fin 
3Z  pr  Figure  S 
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Best  Interpretation  of  f  versus  N^e  of  Tiree  Matrices 
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Schematic  Diagram  of  Apparatus 
Figure  3 
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Skewed  pin  Matrix 
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Recorded  Data 


“  70.0oF.  (a  led  .nt  air  temperature) 


^orifice  -  1+.6 2  in.  HpO 
APnatrix  *  10.2  in.  H2O 
Fa  s  30.21  in.  Hga  (atmospheric  pressure) 
d  =  2,310  in.  (orifice  diameter) 

D  -  3.08  in.  (pipe  diameter) 

ti“tfi  *»  2.73  niv  (initial  temperature  difference,  3  thermocouples) 

chart  speed  s  2,0  sec/in 

chart  scale  s  0.553  mv/ir 

max.  slope  -  -6.772  in/ in 

'ats  *  Oo33U  lhm  (matrix  mass)  (Table  VIII) 

cs  =  0.216  BTU/lbm  °F,  (matrix  specific  heat)  (Table  VIII) 

Djj  *  0,00.35  ft.  (matrix  hydraulic  diameter)  (Table  VIII) 

Aq  ■  0.059li  ft^  (matrix  flow  area)  (Table  VI  Cl) 
a  *  ii3.L0  ft^  (matri<  heat  transfer  area)  (Table  VIII) 

Flow  hate  Determinate  on 


Using  the  method  in  (3),  the  fluid  flow  rate  will  be  found 

=  3 59  K  d“  Fa  Y  orifice,  lbra/hr  |^(3)s  p.  57,  equation  3 

h  =  flow  coefficient,  velocity  of  approach  factor  included 
( dimens ionless) 


Fa  =  thermal  expansion  factor  (dimensionless) 
I  =  expansion  factor  (dimensionless) 


hi 


*  density  of  flu j-c!  flc.:!:  at  Ini' t  :iio  or  orifice  ( 1  ^'P,/ ft' ) 
rtf)  -  nuynold  s  number  based  on  pip?  di.meter 
/&  -  d/D  =  c.75 
assuming  iij-j  ®  9it,000 

K  -  0.7)473,  [(3),  p.  26,  Table  5] 

d^  =  5.33c  in^ 

I  -  0.9959  [(3),  p.  71,  Fig.  UOb] 


/O  =  n  7c;a  t  /ff3  fe  -  -  ons'b  llofy<z 

(  Oo756  lhii/ft-3  ~  (T3^)(FJZr  K 


Fa  -  1.000  j(3),  r.  67,  fig.  33] 

;/f  a  (359) (o.7li73) (5. 336) (1.000) (0.9959)  /  (0,0756)(I*.62) 

*  839  lbm/hr 

Rq  “  0 . 00,1 2lih  Wf/D/U  jO)  p.  58,  equation  o] 

yW  -  1.225  x  10-5  ibm/ft-sec,  at  70  °F,  Rio),  0.  69 
Table  2-8]  L 

Rd  -  ( 0 . 00)+ 2li  1; ) (339 )/( 3.03) (1.225  x  10"?)  «  9)4,326 

The  value  of  *  9U,326  does  not  change  the  initial  flow  coefficient, 
u"f  J  Do  te  r.n  i  na  t  i  o 

To  obtain  NTU,  it  is  necessary  to  determine  the  experi~ 
mental  maximum  slope  of  the  g<  n  ralized  slope  of  the  generalized 
heating  curves 


<%fl  -  tx 


Jt(r/NTC) 


MAH 


i’ron  the  cHfiniti  3 


'Y  =  Jih _ 

\TS  cs 


Q  - 


Wl 

ur? 


and 


NT  U  ^ 


h  A 

cvfCp  ) 


then 


1 7ntu  = 


ft  (0- 


CJ£  Cf 

w 


Wt 


c *Sl 


and 


d{iyNTu)= 


Also 


£  f  1  -  ti 


*fg- .  1 


and 


Thus  the  expression  for  thr.  .naxiunni  slope  b- corses 


pfefr) 

~  fWs  c5"  \  - 1 - 

c$(t f  ?  -  t  J-l) 

dY^HTU) 

v  c  5.  /  £  j.  j  -■£  h 

MAX 

<£  e 

MAK 


./here  tne  quantity 

cflfc  f?  -  t  fi ) 

d  e  . 


MAX 


hi 


ia  the  maxim.- i  lo,  o;  t  ,«<!  .-outlet  luid  trmp  „i  itwi"  iif 
f.rence  versus  tin  ,*v  »  ,.c  tills  t  'iptruturo  dif f-.r^nct  is 

recorded  as  a  fuictio  >i  tins-:,  th  i-.admum  slops  is  *:I.'t«i:n in¬ 
directly  from  th.-.  i  >;p  r:  mental  traces.  3in"-j  the  slops  is 
ineas ur-d  directly  fro  the  chart,  it  *hc  cotr.pu tnttono  it  is 
necessary  to  introduce  the  cii  n’t  scale-factors. 

cf  =  0.210b  BTU/lbn  °F,  fluid  specific  heat  at  JO  °F, 
[(10),  p.  Ijl,  Table  2-3] 

~W l  Cs  _  (O.  3  3  4){o.^lG)(3(s>Oo)/( &'Sq)(o. 24-04) 

ZJT^  ~ 

-  /.  2-  6  9  sec 


(MX^NTU) 


(/.  Z&9) 

-(3.72) 


(-  6.  773 


NTU  *  9.16$  (Table  IX  -  Appendix  II) 


Average  flui d  temperature  determination 

tn  =  70  °F»  ■  1.07  mv 

ip -if 2  “  2,73  mv  (5  thermocouples) 

^i'-tfi  a  1'H  "  2 ‘73  =  0.516  mv  (l  thermocouple) 

5  ~tj"~ 

tx  » (ti«tf ]J+t.f i  “  0.516  +  1,07  B  1.616  mv  s  88.87  hi- 

^ ivg.  *»  ’yi!  71  a  70  *  88.87  *  97,11  55  p 

2  T 

Heat  transfer  Parameter  Determination 

Nst  *  NTb  Ac  *  (9.165) (0.0591)  =  0.0126 

a  m  j(o 

l‘Jpr2/7  »  C.  at  tay  ,  [(10),  p.  Table  2-10  j 

:jdtN?r2/3  =  (0.0126)  (0.795)  =  0.00?;  8 


0.8934- 


rteynolds  Numbjr  Detr.rr:.  i  ati.  i 


5 


[J  s*  U 

77f 


;;-r 


A  r 


(339)(36C'0;/(o.05?2l)  '  3.92  lbn./cac-ft2 


/>i  -  1*223  x  10-2  Ibm/ft-sec,  at  t„w 
Table  2-3] 

•V,.  *  (2.00135) (3.9-' )  /  (1.225  :c  30<)  -  332 


[(10;,  p.  6°, 


i? fiction  factor  De  term  nation 


Iht  friction  data  wars  taken  at  anbi«nt  air  t  -mper.-lar  ss, 

P f  *  -  =  pavg 

vf  R  Tfl 

PavK  *  ‘1  -  ^  matrix 

- 2” 

pavg  -  (30.;  \)(13.6)  -  10.2/2  =  305.9  in.  H20 

Pt  -  (U05.9)(5.20U)/(53.3)(530)  -  .073?  lbm/ft.3 

f  -  (10.2) (5. 203) (0.0737) (2) (32. 2) (0.059i  V  /  (3. 92)2(33. 30) 
=  0,0225 
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SOUJTi'7"  IQ..  T.  "  :Lim!  CLJPZ  OF  THS  .USNiiiULIiCL  uFiffiNG  CURVE 


Locke  (l)  derived  the  expression  for  thr  slope  of  the 
generalised  heating  curve  equal  o:  page  7)>  noting  that  when 

X  -  L  tuat  Z  H  NTU, 


M 


;  7  /'  , - \  .  - 1 

i  1  .  b-v.  /N  T  'Y  /  L 


1  /  * 


l  ' 


a.  JL*C  J-60li  digital  cxnputer  was  usw.  to  jolve  Vm  equation# 
j~  ~  ^  r  •  \  \ 

Values  of  I  y.I  r  were  obtained  by 

•  S'  ~  /  SO'"  •"  .  j 
i  v(  f  /  ^  i 

varying  r/NTU  for  particular  values  of  NTH.  r'or  values  of  NTU  - 
and  using  an  increment  of  !  /NTU  «  ,00001*  the  naximum  slop-  was 
obt.v ntd.  at  the  initial  value  c.f  / VTU  effectively  giving 
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/  M 


K;  T 


J  /'* 
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M  h  * 
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\ 

! 

j 


-or  NTJ  >  2,  1  ncreiurnts  of  ;  /NTU  *  ,0C1  w,r  ••  used 


r /n?it 


increased  until  a  value  of  the  slope  smaller  than  the  previous 
value  was  obtained,  and  this  raxi/mr-  slope  was  recorded.  The  com 


outer  then  moved  to  the  next  value  of  MTU, 


!  1  — 

The  values  of i  '  ' 


■f-)! 


for  values  of  NTU 


/  s 

and  i  -  /NTU  are  cor  J\n  . 


I-  *.  V 


Tabl «[X 


Values  of  the  Max.  Slope  of  tlie  Generalized  Heating  Curve  »NTU  and  "VOTU 
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